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As the main excitation source in the process of vehicle driving, road roughness directly 
affects the vehicle ride comfort and handling stability. In this paper, an on-line estimation 
method of road profile based on vehicle dynamic response is proposed, which estimates 
the road disturbance using the transmission characteristics of the system and realizes 
the road roughness classification by calculating the road information via Power Spectral 
Density (PSD) estimation. Since the existing standard calculation method of PSD of road 
disturbance can only be used to process continuous data, an approximate method suitable 
for discrete data is proposed. In addition, in order to coordinate the response speed and 
estimation accuracy of the system, the mechanism of the influence of time-window on 
the estimation accuracy is studied, and the optimal length of time-window is determined. 
The estimation accuracy under the condition of sprung mass and damping coefficient vari¬ 
ation is also analyzed. Finally, al/4 semi-active suspension test bench is conducted, and 
the results show that the system parameters (sprung mass and damping coefficient) and 
driving speed have little effect on the accuracy of road estimation. At the same time, the 
accuracy of on-line estimation method of road profile under five different levels of road 
is more than 90%, and the overall estimation accuracy rate reaches 98.5%, which verifies 
the effectiveness of on-line road estimation method. 

© 2019 Elsevier Ltd. All rights reserved. 


1. Introduction 

As the main excitation source of the vehicle, road profile directly affects the dynamic performance of vehicles (including 
ride comfort and handling stability), that is, vehicle dynamic control requires to always be based on the driving road. 
Therefore, real-time perception of road roughness is of great significance to improve the performance of control system, 
vehicle ride comfort and handling stability. 

At present, scholars have done some research on road profile estimation. The existing road estimation methods can be 
divided into the following four categories: direct measurement, non-contact measurement, neural network based estimation 
and system dynamic response based estimation. 

The direct measurement method is to measure the road elevation by means of contact with the ground through special 
vehicles [1-7] equipped with road roughness meters (such as profilometer, longitudinal profile analyzer, etc.). Recently, 
Nueraihemaitijiang Abulizi et al. [1 ] designed a road roughness estimator composed of a compact road profiler and ArcGIS, 
which is used for road management instead of vehicle system. Nikolay Lushnikov and Petr Lushnikov [2] compared the 
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characteristics of several profilometers used for measuring of longitudinal roughness of the road surface and studied the 
accuracy of road roughness measurement with them. Because of the poor transplant ability, this method is mainly used 
for the maintenance of transportation roads, and cannot be applied to real-time measurement of conventional vehicles. 
The noncontact measurement method obtains the road elevation by scanning the road surface using special sensors installed 
on the vehicle body, such as laser sensors and ultrasonic sensors [8-10], Hyo-Jun Kim et al. [8] designed a road profile mea¬ 
suring system based on the computational road-sensing algorithm which adopted a composite-sensor system with two laser 
sensors. Yong Li et al. [9] presented a road detection method using airborne LiDAR point clouds which obtained a correctness 
of 76 percent confirmed by practical data. Compared with the direct measurement method, this method has better trans¬ 
plant ability, and can be used for real-time detection of general vehicles. However, the cost of special sensor is expensive, 
and its measurement effect is limited by the environment [11,12], In the rain and snow weather, the accuracy of measure¬ 
ment is difficult to be guaranteed. 

In the neural network based estimation method, Ali Solhmirzaei et al. [13] constructed a multi-input and multi-output 
feed forward wavelet neural network trained by the BP algorithm to estimate the road profile. Similarly, H.M. Ngwangwa 
et al. [14,15] and Mahdi Yousefzadehal et al. [16,17] realized the reconstruction of road defects and road roughness based 
on the neural network which trained by the dynamic response signals (acceleration of sprung mass) of vehicle. Yechen 
Qin et al. [18] proposed a road estimation method with application of adaptive Neuro-Fuzzy Inference System (ANFIS), which 
integrates an improved distance evaluation technique and a 2-layers ANFIS classifier. All the methods above only studied the 
measurement of road profile with statistical characteristics, but road profile also has instantaneous characteristics. H.M. 
Ngwangwa et al. [19] then presented a road profile estimation method which performed particularly well in the estimation 
of road damage (such as a bump, depression, or pit, etc.). The test results based on the heavy truck and small underground 
mine truck showed that the estimation error of this method is less than 25%. Based on the advantages of neural network 
technology, the problems of transplant ability and equipment cost are solved and a solution for the development of 
vehicle-borne high-performance road estimation method is provided. However, the creation and training of artificial neural 
network require a large number of original data which requires to be collected under different roads in advance. In addition, 
in the process of road estimation, this method occupies a large amount of computing resources and a long response time of 
the system. 

Similar to the neural network based method, the system dynamic response based method does not need external special 
measuring device, and has the advantages of good transplant ability and low cost. The difference is that the system dynamic 
response based method dose not depend on prior knowledge and has less computational load, so it has a broader research 
prospects. P.E. Uys et al. [20] developed a road profile identification method based on the Fourier analysis of displacement of 
Power Spectral Densities. Fei Hu et al. [21 ] proposed a fast, computationally efficient, and economical road surface recovery 
method based on a vehicle dynamic model with output tracking feedback control, which only adopting accelerometers. Sim¬ 
ilarly, Keum-Shik Hong et al. [22] presented an estimation method of road condition using the body acceleration based on 
the developed linearized new model and according to which, a modified skyhook controller is proposed. Moustapha Doumi- 
ati et al. [23] used the measured vertical acceleration to calculate the vehicle body position and based on which, the road 
profile and the tire load are estimated with the Kalman filter. Yuchuan Du et al. [24] developed an International Roughness 
Index (IRI) measurement system using a Z-Axis acceleration and a GPS device, and the measurement accuracy is verified by 
single-wheel linear model and two-wheel linear model. M Burger et al. [25] proposed a road profile and road roughness com¬ 
putational method based on the quarter-car model and half-car model, which depends on the vertical acceleration mea¬ 
sured. Boyu Zhao et al. [26] developed a road profile estimation algorithm by measuring the dynamic response of the 
vehicle with a smartphone. Benefits from the Kalman filter and Rauch-Tung-Streiber (RTS) smoothing technology, the sat¬ 
isfactory estimation accuracy is implemented. Although the system dynamic response based method effectively solves 
the above problems, it still has several shortcomings, such as the influence of the variation of parameters on the estimation 
effect has not been considered. But T. Heyns et al. [27 ] proposed a method for on-line road profile identification of haul roads 
based on Bayesian parameter estimation, which used in road management instead of vehicle system. And a speed calibration 
function is developed in this method to eliminate the influence of speed. Similarly, Yechen Qin et al. [28] presented a speed 
independent road level classification method, which can effectively estimate the current road level through on-line and off¬ 
line estimation. Wuhui Yu et al. [29] utilized a linear Kalman filter to online estimate the road roughness, and considered the 
influence of load variation, the recursive least-squares estimation method was adopted to calculate the current sprung mass 
and based on which, the estimation accuracy was effectively improved. Yechen Qin et al. [30] then applied the pattern recog¬ 
nition method commonly used in fault diagnosis to road profile estimation and designed a random forest classifier based 
road profile estimator, which taken the damping variation into consideration. In the above research, scholars began to realize 
the impact of the variation of parameters on the estimation accuracy, and studied the impact of a single parameter and based 
on which, the stability of the road estimation system was improved. However, all the above studies are based on the passive 
suspension, with the development of semi-active and active suspension, many parameters of system may change. There are 
few studies relate to the road estimation technology under the condition of time-varying parameters. 

According to the existing shortcomings, a novel on-line estimation method of road profile based on unsprung mass accel¬ 
eration signal is proposed in this paper, which mainly includes road disturbance estimation and road roughness classifica¬ 
tion. In the road disturbance estimation, a road disturbance estimator based on the transmission characteristics of 
suspension system is designed. According to the measured road disturbance, a road roughness classification method based 
on the Power Spectral Density (PSD) estimation is proposed. Furthermore, the mechanism of the influence of time-window 
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on the estimation accuracy is studied, and the estimation accuracy under the condition of sprung mass and damping coef- 
ficient variation is analyzed. Finally, a 1/4 semi-active suspension test bench is conducted, and the effectiveness of this 
method under the condition of system parameters (i.e. sprung mass and damping coefficient) and driving speed variation 
is verified. 

The rest of this paper is organized as follows: Section 2 selects the 1/4 semi-active suspension as the research object and 
establishes its numerical model. Section 3 designs a road disturbance estimator based on the unsprung mass acceleration 
signal, and proposes a road roughness classification method based on the PSD estimation. Section 4 studies the mechanism 
of the influence of time-window on the estimation accuracy, and analyzes the estimation accuracy under the condition of 
parameters (sprung mass and damping coefficient) variation. Section 5 conducts a test bench and verifies the effectiveness 
of the proposed method. Section 6 presents the general conclusion. 

2. System numerical model 

The on-line road profile estimation method proposed in this paper is a model-based method and the dynamic response 
characteristic of the vehicle is the research basis of this method. To make the method simple and easy to implement, a 1 /4 
vehicle linear model is selected as the research object. Although this model is simple, it captures the most basic feature of the 
vertical model of the vehicle. The system model is shown in Fig. 1 . 

According to the Newton’s second law, the dynamic differential equation of the semi-active suspension is established. 



The road input is filtered white noise, which obeys the Gauss orthogonal distribution. 



where m s and m t are the sprung mass and unsprung mass, respectively: c is the variable damping coefficient: k s is the spring 
stiffness: k t is the tire stiffness; /„ is cut-off frequency: v is vehicle speed; w(t) is time domain signal of white noise; G,(n 0 ) is 
the road roughness coefficient (the specific values are shown in Table 1). 

The control objectives of vehicle under different road conditions are different [31 ]. As shown in Fig. 2 (taking the B, C and 
D classes of roads as an example), when the vehicle travels on the flat road (B-class road) with high speed, it is necessary to 
ensure that the wheels have good contact with the road surface, so as to improve the handling stability. But when in bumpy 
road (D-class road), the driving speed is often low. At this time, ride comfort is the main control objective, and it is significant 
to ensure that the suspension working space within a reasonable range, that is, the wheel will not impact the limit block. 
Because the requirements of vehicle dynamic performance on different levels of road are different, the suspension control 
base on road information is of great significance, which is also the starting point of this paper. 

3. Design of on-line estimation method of road profile 

As can be seen from the previous section, the control objectives of the vehicle dynamic performance depend on the cur¬ 
rent level of road. Therefore, taking the semi-active suspension as the research object, an on-line road estimation method 



Fig. 1. Simplified 2-DOF 1/4 vehicle linear model. 
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Table 1 

Parameters for a 1/4 vehicle model. 

Parameter Value 


Sprung mass m s 
Unsprung mass m t 
Damping coefficient 
Spring stiffness k s 
Tire stiffness le t 


317.5-635 kg 
45.4 kg 

1000-7000 N-s/m 
20000 N/m 
192000 N/m 


Road Excitation 


Unsprung Acceleration 



Road Disturbance 


Power Spectrum 
Density 
Calculation 


Variable 


1 

1 

Interval 


1 

1 

Output 

parameters 

Cmr&M„ 1 
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Judgment 

V_ 


! 

Road Level 


Fig. 2. Structure of the proposed estimation method. 


including road disturbance estimation and road roughness classification is designed, and its overall structure is shown in 
Fig. 2. Firstly, an estimator of road disturbance is designed and then based on which, a road roughness classifier is developed. 


3.1. Estimation of road disturbance 


In this paper, a road disturbance estimator is designed based on the transmission characteristics of the suspension sys¬ 
tem, and the sprung mass acceleration is selected as the measurement object. 

X(s)-m (3) 


By making Laplace transformation of formula (1) and eliminating the z s (s), we can obtain: 

( m t s 2 +cs + k s + kt)Zt(s) — (fts + cs) ^ = kz ^ 

v ' w m s s 2 + cs + k s w w 


(4) 


By making Laplace transformation of formula (3) and substituting it to the formula (4), the transfer relationship between 
wheel vertical acceleration and road input can be obtained. 

Zri^l _ b 4 s 4 + b 3 s 3 + b 2 s 2 + b^s + b 0 /si 

x(s) = Q4S 4 + a 3 s 3 + a 2 s 2 + cls + a 0 


where a 4 = m s m t , a 3 = k t c, a 2 = k s k t , a, - a„ - 0, b 4 = m s m t , b 3 = m s c + m t c, b 2 = m s k s + m t k s + m s k t , b, = k t c, b 0 = k s k t . 

Then, the numerical model of the system, which takes the acceleration signal as input and road disturbance as output, can 
be expressed as follows: 

a 4 z' 4) (t) + a 3 z^(t) + a 2 z r (t ) = b 4 x (4) (t ) + b 3 x (3) (t ) + b 2 x(t ) + b 4 x(t) + b 0 x(t ) (6) 

From the expression of the transfer function of system, it can be seen that the road disturbance can directly cause the 
vertical vibration of wheels. Thus, the road elevation can be obtained indirectly by detecting the unsprung mass acceleration 
during the vehicle driving process with acceleration sensor. 

Based on the above formulas, a simulation model based on the system transfer function is established, and the system 
parameters are set according to Table 1. Fig. 3 shows the disturbance of C-class road generated by formula (2) when the 
driving speed is 20 m/s. Fig. 4 gives the measured vertical acceleration signal of unsprung mass in the time domain. Fig. 5 
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Fig. 3. Input road disturbance. 



Time/s 


Fig. 4. Measured unsprung mass acceleration. 


shows the comparison between the road disturbance estimated by the dynamic response model and the original input dis¬ 
turbance in the time domain. 

It can be seen from the Fig. 5 that the road information obtained based on the dynamic transfer characteristics of the sys¬ 
tem can well track the original input road, which verifies the validity of the designed road disturbance estimator and the 
accuracy of the established system model, and lays a foundation for the road classification later. 

3.2. Classification of road roughness 

In the previous part, the road disturbance is obtained based on the transmission characteristics of the system. However, 
the level of road roughness is ultimately required for the suspension control. At present, PSD is commonly adopted to 
describe the statistical characteristics of road roughness [32], Thus, the analysis of PSD of road disturbance is necessary 
to be carried out. A large number of experiments show that the vertical acceleration data belong to a stationary random pro¬ 
cess with trend term. Therefore, the AR model with good variance performance and high spectral resolution is more suitable 
for processing this type of data [33], 

The key to the establishment of AR model is the solution of parameters. In order to reduce the computational difficulty 
and shorten the computational time, the Burg algorithm [34-36] is selected. Burg algorithm skillfully avoids the calculation 
of autocorrelation function and can realize the estimation with fewer data records. Furthermore, the estimated results based 
on this method are very close to the real values. 

Assuming that the road roughness z r (k) in the frequency domain satisfies the formula of AR model: 

Zr(fc) = -CjZ r (i - 1) - c 2 z r (i -2) - c p z r {i - p) + e„ (7) 

where c p (p = 1,2, ■ • •, p) is the model parameter; e d is the stationary white noise with zero mean and a\, variance. 

Assuming that the forward and backward prediction errors of linear predictive AR model ar ef p (k) and b p (k), respectively, 
f f„(k) = Zr{k) + a p iZ r {k - 1) + ■ ■ • + a pp z r {k - p) 

\ b p {k) = z r (k -p) + a„iz r (k-p + !) + ••■ + a pp z r {k ) 
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Fig. 5. Comparison between estimated and input road disturbances. 


Then, the sum of the power of the forward and backward prediction errors is: 


p fl> - 2 IP/ + P&] - jvTTj, S 11 2 + 5Z M k ) I 2 

/ p (k) and fa p (k) have the following recursive relationships: 

f/ s (fe)=/ s -,(fe) + W s -i(k-l) 

\b s {k) = fas-i(k) + h_sf s _^(k - 1) 


(9) 


(10) 


where s is the order, s = 1,2, ■, p; h s is the reflection coefficient. 

According to the Burg method, and making the sum of energy of forward and backward prediction error the smallest rel¬ 
ative to the reflection coefficient, the parameters of AR model can be obtained based on the Levinson-Durbin recursive algo¬ 
rithm [37,38], 

Fig. 6 shows the comparison between the PSD of the estimated road disturbance and that of the input road disturbance, in 
the spatial frequency domain. As can be seen from the figures, the analysis results in the frequency domain are basically con¬ 
sistent with those in the time domain, which verifies the correctness of the results in the time domain and shows the fea¬ 
sibility of adopting this method to classify the road roughness. 

Based on the PSD of road disturbance, the road roughness requires to classified. The spatial frequency of common road 
ranges from 0.011 m J to 2.83 m ', according to the ISO 8601 [39], the PSD of the standard road can be fitted by the following 
formula: 


G,(n)-G,(n 0 )(£) “ (11) 

where n is the spatial frequency in m 1 ; n 0 is the reference frequency in m 1 ; w is the frequency exponent which usually 
equals to 2; G q (rio) is the road roughness coefficient, and its value depends on the level of road. 



Spatial frequency/m" 1 


Fig. 6. PSD comparation between estimated and input road disturbances. 
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Based on the formula (11), the Root Mean Square (RMS) value of PSD of road disturbance for each level of road in a range 
of 0.01 lm 1 to2.83m 1 can be obtained as shown in Table 3. It can be seen that the RMS value of PSD of the same type of road 
is distributed in a specific range. Therefore, by calculating the RMS value of PSD of estimated road disturbance and judging its 
interval, the level of the road roughness can be detected. It should be noted that the formula (11) is only suitable for the 
calculation of continuous data, but the signal measured by sensors are discrete. Thus, an approximated calculation method 
based on the discrete data is proposed as shown in formula (12): 

a * = Y< + G i( n ( x + 1 >))(»(*) - "(* + !)) O 2 ) 

where x is the number of data points collected by the designed road estimation system in the time-window t; i and j are data 
point numbers corresponding to spatial frequencies n = 0.011 and n = 2.83, respectively: n(x) is the spatial frequency corre¬ 
sponding to data number x; G q (n(x)) is the PSD of road disturbance when the spatial frequency is n(x). 

Based on the formula (12), the RMS value of PSD of road disturbance in current time-window can be calculated, and then 
the level of the current road can be detected by judging the interval of the value according to the Table 3. 


4. Sensitivity analysis 

4.1. Relationship between estimation accuracy and time-window 

From the calculation method of PSD of road disturbance proposed in formula (12), it can be seen that the estimation accu¬ 
racy of on-line estimation method of road profile is directly related to the accuracy of calculation method. However, the 
accuracy of the calculation method is affected by the length of the time-window during the calculation process. Therefore, 
in order to obtain higher estimation accuracy, the potential relationship between estimation accuracy and time-window 
requires to be explored. Under the conditions set in section 3, and setting time-window t to be 1 s, 2 s, 4 s, and 8 s respec¬ 
tively, the change trend of RMS value of PSD with time-window is obtained as shown in Fig. 7. 

The “Upper limit” and “Lower limit” in the figure are the upper and lower limits of the RMS value of PSD of C-class road in 
Table 2. According to the road roughness classification method mentioned above, if the RMS values of PSD calculated based 
on the time-window t fall between the “Upper limit" and “Lower limit”, the road passing in this time period can be judged to 
be C-class road. It can be seen from the figure that the RMS values of PSD fluctuate more intensely in a narrower time- 
window, because the uncertainty of the PSD values under each frequency corresponding to the data covered by time- 
window is larger [40], But with the extension of the time-window, the RMS values tend to be stable. When the length of 
time-window is 1 s, many results fall outside the interval of C-class road, i.e. judgment errors occur in the system. With 
the increase of time-window length, the RMS values of PSD begin to converge between the upper and lower limits i.e. the 
estimation accuracy is gradually improved. Furthermore, it can be seen that the extension of time-window obviously 
increasing the response hysteresis of RMS values on road roughness, which presents as a curve with a smaller slope at 
the beginning of the figure. 

From the above research, it is significant to select a reasonable length of time-window when adopting the RMS value of 
PSD of road disturbance to classify the road roughness. But there is a contradiction in the selection of time window, that is, 
the contradiction between the estimation accuracy and response time of the system: Extending the time-window can 
improve the estimation accuracy, but it will increase the response time, which will make the system unable to complete 



1000 1500 

Distance along the road/m 


Fig. 7. RMS values of PSD of road disturbance under different time-window lengths. 
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Table 2 

Control objectives of different road conditions. 

Road roughness level Road roughness Driving speed Control objective 

B 64xl(T 5 100 km/h Road holding 

C 256 x 10 5 60 km/h Compromise 

D 1024 x 10 -5 30 km/h Ride comfort 


Table 3 

Classification standards of road roughness. 

Road roughness level Value G,(n 0 )/m 3 

A 8 — 32 x 10 -6 

B 32 — 128 x 1(T 6 

C 128 — 512 x 10~ 6 

D 512 ~ 2048 x 10 6 

E 2048 ~ 8192 x lO 6 

F 8192 ~32768 x 10 6 

G 32768 ~131072 x ltT® 

H 131072 ~524288 x1(T 6 


2.69-5.38 x 1<T 3 
5.38 - 10.77 x 10 3 
10.77 — 21.53 x 10 -3 
21.53 ~ 43.06 xl0“ 3 
43.06 ~ 86.13 x 10“ 3 
86.13 - 172.26 x 10“ 3 
172.26 — 344.52 x 10~' 
344.52 - 689.04 x 1CT : 


the judgment of road level in a relatively short time; Shortening the time-window can improve the response speed, but it 
will deteriorate the estimation accuracy and deteriorate the effect of subsequent vehicle control. Thus, in order to determine 
the optimal length of time-window, the estimation error rate of the system under the corresponding time-window is 
selected as the evaluation index of estimation accuracy. 

E r = — x 100% (13) 

Pt 

where the E r is the estimation error rate of the system; p t is the number of total estimated points; p e is the number of incor¬ 
rectly estimated points. 

In the actual driving process, since the A ~ E classes of roads are the most common and important and F ~ H classes of 
roads are rarely to see, this paper mainly focuses on the analysis of these five roads. Selecting time-windows from 1 to 
10 s in step length 0.1 s, the estimation error rate of A ~ E classes of roads under different time-windows can be obtained 
as shown in Fig. 8. 

From the Fig. 8, it can be seen that no matter what level of road the vehicle is driving on, the estimation error rate of the 
system decreases obviously with the increase of the length of time-window i.e. the estimation accuracy significantly 
improves. But when the time-window increases to a certain length, the change of the estimation error rate tends to be gentle, 
that is, extending the time-window is of little significance to improve the estimation accuracy. On the contrary, it will 
increase the response time and deteriorate the response speed of the system. In the figure, the gray marks the points with 
estimation error rate of 0% under each level of road, and the corresponding lateral axis numbers are the minimum time- 
window lengths required by the on-line estimation method of road profile to distinguish each section. In operational appli¬ 
cation, the time-window can be chosen as long as 4.8 s, 5.2 s, 5.6 s, 6.3 s, and 8.4 s for accurate classification according to 
Fig. 7. The red “°” marks the point with the variation rate of error rate approaching 0% under E-class road. It can be seen that 
although the error rate of this point is not 0%, the error rate has hardly changed since that point when the time-window con¬ 
tinuous to be extended. In order to balance the accuracy demands of control systems and tolerance of delay, the 7.3 s of time- 
window length is selected under E-class road, which is also selected as the final time-window length in this paper to discover 
the accuracy performance of the method. 

4.2. Analysis of the variation of RMS values of PSD with suspension parameters 

In the above contents, the calculation method of PSD of road disturbance is proposed, the rules of estimation accuracy and 
respond speed varying with time-window are studied, and the selection method of time-window length is clarified. But it 
can be seen from the formulas (5) and (12) that the calculation accuracy of PSD of road disturbance is not only related to 
the input road and time-window, but also to the transmission characteristics of the system. The transmission characteristics 
of the system depend on its own parameters, including sprung mass, unsprung mass, spring stiffness and damping 
coefficient etc. 

With the popularization of semi-active suspension, only the sprung mass, spring stiffness and damping coefficient 
can be changed in vehicle suspension systems. Among these three parameters, the variation of spring stiffness can only occur 
in the semi-active suspensions with air spring, while the variable damping coefficient is easy to achieve and common in 
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Fig. 8. The results of minimum time window selection. 


semi-active suspensions. In addition, because of the variable number of passengers, the variation of sprung mass is very 
large, which is also significant. Therefore, this paper only studies the relationship of RMS values of PSD with sprung mass 
and damping coefficient, the influence of variable spring stiffness in semi-active suspension with air spring on estimation 
accuracy is temporarily overlooked. The RMS values of PSD of road disturbance under A ~ E classes of roads are calculated 
in the range of 317.5 kg ~ 635 kg of sprung mass. This range is widely enough to cover all the possible sprung mass of vehicle. 
The calculation results are shown in Fig. 9. 

Where M 0 is initial sprung mass 317.5 kg; M is current spring mass; RMS 0 is the root mean square value of PSD of esti¬ 
mated road disturbance under initial sprung mass; RMS is the root mean square value of PSD of estimated road disturbance 
under current sprung mass. 

From Fig. 9, it can be seen that the variation of sprung mass has a little influence on the RMS values of PSD. Although the 
sprung mass increases to the double of the initial value, the error between RMS value under current sprung mass and that 
under initial sprung mass is less than 1%. Because the influence of sprung mass on RMS values of PSD is small, the system can 
adapt to the variation of sprung mass, further research is not carried out on this topic. 

Considering that the damping coefficient varies greatly with road condition, the range of damping coefficient is set to be 
1000-7000 N ■ s/m. Within this range, the RMS values of PSD of estimated road disturbance are calculated under A~E 
classes of roads. The relationship between RMS values of PSD and damping coefficient is shown in Fig. 10. 

Where Co is initial damping coefficient of suspension 1000 N ■ s/m; C is current damping coefficient. 

It can be seen from Fig. 10 that, with the increase of damping coefficient, the RMS values of PSD estimated by system 
continuously decrease. Furthermore, with the improvement of road level i.e. the increase of road elevation, the influence 
of damping coefficient on RMS values of PSD gradually weakens. The estimation error of system under A-class road is the 
greatest. When damping coefficient varies in the given range under this road, the maximum error between RMS values of 
PSD under current damping coefficient and that under initial damping coefficient is only 0.73%. According to the simulation 
results, it can be seen that the variation of damping coefficient also has little effect on estimation accuracy of road estimation 
system, which can be neglected. 



M/M q /kg 


Fig. 9. How sprung mass affects RMS values. 
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C/Cg / N-s/m 

Fig. 10. How damping coefficient affects RMS values. 


In this paper, the influence of two key variable parameters of semi-active suspension on estimation accuracy of the sys¬ 
tem is studied. The simulation results show that no matter which parameter changes, the maximum error between RMS val¬ 
ues of PSD estimated under this condition and initial value is less than 1%, that is, the variation of parameters in semi-active 
suspension has little effect on the estimation accuracy of the designed road estimation system. The system has good robust¬ 
ness and can satisfy the requirements of semi-active suspension. 

5. Bench test 

In order to verify the effectiveness and stability of the on-line estimation method of road profile under the condition of 
variable parameters, a test platform for 1/4 semi-active suspension is built. 

5.1. Handware-in-the-loop simulation platform 

Fig. 11 shows the 1/4 semi-active suspension test bench. MR (Magnetorheological) damper with variable damping coef¬ 
ficient is adopted as the actuator in this test bench, and the sprung mass is adjusted by increasing or decreasing the number 
of sandbags. The whole test bench includes host computer, controller, MR damper driver, INSTRON 8800 numerical control 
hydraulic servo exciter, acceleration sensor and data collector. The functions of each module are as follows: 

(1) The host computer is used to establish the simulation model in Matlab/Simlink software, set simulation conditions, 
observe the waveforms of various physical signals in real time and analyze and save data of sensor; 

(2) The controller is the rapid prototype controller based on D2P (Development to Prototype), which is used for Running 
control algorithm and output control voltage signal to MR damper driver; 



Fig. 11. Configuration of handware-in-the-loop simulations. 
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(3) The MR damper driver is used to receive the PWM signal from the controller and drive the MR damper, at the same 
time, feeds its control current signal back to the D2P to realize the calculation of inverse model of MR damper, so as to 
complete the closed-loop control; 

(4) INSTRON 8800 numerical control hydraulic servo exciter is used to receive signals and simulate the displacement of 
suspension, and transmit the force signal of MR damper to the rapid prototyping controller; 

(5) The acceleration sensor is used to measure the acceleration signal of unsprung mass at real-time; 

(6) The data collector is used to receive the sensor signal and transmit it to the controller and host system. 

5.2. Signal processing 


According to the section 2, obtaining accurate road elevation signal in time domain is the key to realize the on-line esti¬ 
mation of road profile, while obtaining accurate unsprung mass acceleration signal is the premise of realizing estimation of 
road elevation. Thus, in order to prevent the interference of the high-frequency noise caused by engine vibration during the 
acquisition of the unsprung mass acceleration signal and solve the drift of the accelerometers due to the DC offset, a com¬ 
bined filter is designed in this paper to process the acceleration signal as shown in Fig. 12. The combined filter consists of a 
first-order low-pass filter I(s) and a first-order high-pass filter H(s). Among them, the first-order low-pass filter L(s) is 
adopted to filter the high-frequency noise as well as the first-order high-pass filter H(s) is utilized to eliminate the drift 
of the accelerometers. The transfer functions of the low-pass filter and high-pass filter are as follows: 


L(s) = 


GlCQl 
co L + s 


(14) 


H(S) mi 
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COh + S 

Then, the transfer function of the combined filter can be expressed as follows: 
sw l G l G h 
(co l + s)((Oh + s) 


H(s)=t- 


(15) 


(16) 


where a> L and (Oh are the cut-off frequencies of the low-pass filter and high-pass filter, respectively; G t and C H are the gains 
of the low-pass filter and high-pass filter, separately. The parameters involved in Formulas (14) and (15) are shown in Table 4, 
which is determined by numerous tests. 

Additionally, in the process of measurement and transmission of unsprung mass acceleration signal, abnormal data will 
appear because of interference and accidental errors. These erroneous data mixed with a large amount number of correct 
data will cause great errors in the processing results, so it requires to be eliminated in the process of signal processing. Con¬ 
sidering the large number of signals collected in each time-window and the complexity of the processing method, Letts cri¬ 
terion (3er criterion) is adopted to process the abnormal data, that is, if the measurement data x d satisfies | v d \ = \x d - x\ > 3S, 
it is considered that x d is abnormal data and can be removed, otherwise, it will be retained. S is the standard deviation cal¬ 
culated by Bessel formula, and its expression is shown formula (17): 
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Fig. 12. Signal processing. 


Table 4 

Parameters of combined filter. 

Parameter Value Parameter Value 

G l 2 co L 50 

C„ 2 (Oh 0.5 
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where S is the standard deviation of sample data; xfi = 1,2, • • n) is the sample data; x is the average value of sample data; n 
is the number of sample data. 

5.3. Results analysis 

In order to validate the performance of the on-line estimation method of road profile, two groups of tests are conducted in 
this paper. Firstly, the estimation effect of road disturbance estimation method based on the vehicle dynamic response on 
road excitation is studied; Then, the classification results of road roughness are given based on the analysis of estimated PSD. 

5.3. f. Performance of on-line road profile estimator of road disturbance 

In this part, the ability of the on-line road profile estimation method to estimate the road disturbance is test and verified. 
In order to demonstrate the superiority of this method, the algebraic estimation-based estimation method proposed by Mar- 
oua Haddar [41] is chosen for the comparison. 

Test #1; C-class road disturbances, variable damping coefficient 

In this test, C-class road is taken as the input disturbance and the vehicle speed is set as 72 km/h. The driving time is 30 s 
under each condition, and the suspension damping is adjusted by the MR damper. 

Fig. 13(a) shows the variation process of suspension working space with road disturbance. It can be seen that the varia¬ 
tion trend of suspension working space is contrary to that of damping coefficient, and the suspension working space grad¬ 
ually increases with the decrease of damping coefficient. However, the road disturbance estimator presents good robustness 
to the variation of suspension damping. Fig. 13(b) shows the comparison of the estimation results of road disturbance 
between the designed method and the existing road disturbance estimation method. As can be seen in the figure that both 
two methods accurately estimate the road disturbance, and the estimation results are not affected by the variation of sus¬ 
pension damping. But, the estimation accuracy of the road disturbance estimation method designed in this paper is obvi¬ 
ously higher than that of the existing method. To quantitatively analyze the estimation accuracy of the two methods, the 
estimation errors of the two methods are studied and the root mean squared error (RMSE) between the real and estimated 
road is further compared. The comparison results are shown in Table 5. As can be seen from the table, the RMSE of the 
designed road disturbance estimation method is 0.0029, which is significantly less than 0.0115 of the existing method, that 
is, the road disturbance estimation method designed in this paper is superior to the existing method. 

Test #2: C-class road disturbances, different sprung masses 

This test also takes C-class road as the input disturbance, and vehicles keep 72 km/h speed on the road for 90 s, i.e. the 
driving distance is 1800 m. The MR damper works with a fixed damping coefficient (without controller). 

Fig. 14(a) shows the comparison of the estimation results of road disturbance between the designed method and the exist¬ 
ing method under no-load condition (for convenience of comparison, the test results of 75 s ~ 90 s are selected). It can be seen 
that under this condition, the estimation results of the two methods both can well track the input road. However, the estima¬ 
tion accuracy of the existing method is obviously less than that of the method designed in this paper, and the specific estima- 
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Table 5 

Comparison between the designed method and existing method using root mean square error (RMSE) criteria. 


Estimation method 

Root me; 

in squared error RMSE - 


Test #1 

Test #2: No-load Test #2: Full-load 

Designed method 

Existing method 

0.0115 

0.0029 

0.0121 0.0134 

0.0047 0.0064 



Fig. 14. Analysis of performance of on-line road estimator with different sprung mass, (a): Estimation performance of road disturbance (No-load), (b): 
Estimation performance of road disturbance (Full-load), (c): Error of estimation (Full-load). 

tion errors are shown in Table 5. Fig. 1 4(b) gives the comparison of the estimation results under full-load condition, compared 
with that under no-load condition, the estimation accuracy of the two methods both decreases slightly. But, the RMSE of the 
designed method is still less than that of the existing method, which is 0.0064 and 0.0134 separately. The comparison of the 
estimation error of road disturbance under the full-load condition is shown in Fig. 14(c). As can be seen from the figure that 
the RMSE of the designed road disturbance estimation method is about 0.004 m. Although the estimation accuracy is slightly 
worse than that under No-load condition, the degree of deterioration is small enough to satisfy the accuracy requirements of 
subsequent road roughness classification. The estimation error of the existing method is about 0.02 m, which is significantly 
more than that of the designed method, that is, the existing method is inferior to the designed method. 

According to the above test results, the road disturbance estimation method designed in the paper is obviously superior to 
the existing method under the condition of damping and unsprung mass variation, which validates the effectiveness and 
superiority of the designed road disturbance estimation method. 


5.3.2. On-line road classification 

In this part, the tests are mainly carried out to verify the ability of the road roughness classification, and the algebraic 
method based on the time-window developed by Zhongxin Li [40] is selected as the comparison objects. 
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Test #1; Road classification of a random sequence of ISO road profile (at constant vehicle velocity) 

The purpose of this test is to validate the classification performance of the designed on-line estimation method of road 
profile. In order to evaluate the adaptability of on-line estimation method of road profile to the random sequence of road, 
a random sequence covering five ISO road profile disturbances is established in this test. The random sequence has duration 
of 1460 s at a constant vehicle velocity (36 km/h) i.e. representing 14600 m of path; the semi-active suspension works with 
constant damping (without controller). The comparison results are shown in Fig. 14. 

Fig. 15(a) illustrates the estimation results of the random sequence, it can be seen that in the all cases, the estimation 
follows the input experimental road signal. Based on the estimated road disturbance signal in Fig. 15(a), Fig. 15(b) gives 
the corresponding RMS values of PSD of road disturbance calculated by the road roughness classifier. The red dotted line 
from bottom to top in the figure marks the boundary of the RMS value of PSD of A ~ E classes roads given in Table 2. As 
can be seen from the figure that in general, the classification effect of on-line estimation method of road profile is good 
for all roads, only with some misclassifications among the neighbor roads. By judging the interval of RMS values of PSD, 
Fig. 15(c) presents the classification results of test roads obtained by on-line estimation methods, which is compared with 
the classification results of the method described in reference [40], It can be seen from the figure that both two methods can 
classify all road in the given random sequence. But, the classification accuracy of the designed road classification method is 
noticeably higher than that of the existing method, and the probability of correct classification reaches 98.5%. Simultane¬ 
ously, it is not difficult to find that both two methods are prone to produce classification errors under sudden change of road 
conditions, which is due to the large fluctuation of the data in the time-window during the sudden change of road conditions. 
The specific comparison results are shown in Table 6. 

Test #2: Classification of rough road profile at variable vehicle velocity 

In the previous test, the estimation performance of road disturbance and classification performance of road roughness are 
verified at constant vehicle speed. However, during the actual driving process, due to the variation of road conditions, the 



Fig. 15. On-line estimation and classification of random road, (a): Estimation performance of road disturbance, (b): On-line estimation of standardized i 
profile, (c): On-line classification of standardized road profile. 
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Table 6 

Comparison of road classification accuracy between designed method and existing method. 
Estimation method Time-window length Estimation accuracy 

Test#l Test#2 

Designed method 7.3 s 98.5% 96.25% 

Existing method 7.5 s 92.8% 92.21% 


vehicle will not always travel at a constant speed. For example, under the flat road, the driving speed is generally high 
because of the good road condition, while under the rough road, considering the safety of vehicles, the driving speed is 
low. This test is used to evaluate the adaptability of on-line estimation method of road profile to the vehicle speed. At 
the beginning, the vehicle is driven on the A-class road at high speed (72 km/h) and after 146 s, the vehicle passes on D- 
class road at the same speed; finally, at t = 365 s, a step increment of the vehicle speed (from 72 km/h to 36 km/h) is con¬ 
sidered keeping the same road class D. 

Fig. 16(a) shows the test results under the given conditions. It can be seen that under the both two driving speeds, the 
road information estimated by on-line estimation method of road profile can well follow the input road disturbance. How¬ 
ever, because of the sudden change of working conditions, estimation errors appear in the on-line estimation method of road 
profile in a short period when the road profile and vehicle speed are switched. As we can see from the Fig. 16(b) that the RMS 
values of PSD is inside the road class D area independently of the vehicle speed (72 km/h to 36 km/h). Fig. 16(c) shows com¬ 
parison of the road classification results between the designed method and the existing method, and the classification results 
of the designed method present the same classification errors the previous results. However, the convergence rate of the 
road classification results of the designed method is obviously higher than that of the existing method. In the second 
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Fig. 16. On-line road estimation and classification under variable vehicle velocity, (a): Estimation performance of road disturbance, (b): On-line estimation 
of standardized road profile, (c): On-line classification of standardized road profile. 
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time-window after the sudden change of speed, the system re-tracks the input road quickly, and the probability of correct 
classification in this test is 96.25%. Unfortunately, the existing method requires 2 to 3 time-windows to re-track the input 
road, and the probability of correct classification is only 92.21%. The comparison results are shown in Table 6. According 
to the above results, it can be seen that the system can still realize the classification of road roughness under the condition 
of variable speed. 

The above two groups of tests validate the estimation and classification performance of the on-line estimation method of 
road profile, the results show that under the condition of variable environment (including variation of damping coefficient, 
sprung mass and vehicle speed), the designed on-line estimation method of road profile has good adaptability and still can 
accurately realize the estimation and classification of current road. 

6. Conclusion 

In this paper, a novel algorithm of road profile estimation and classification for semi-active suspensions has been pro¬ 
posed. A road disturbance estimator based on the transmission characteristics of suspension system is used to rebuild the 
road profile signal. A road roughness classification algorithm, according to the Power Spectral Density of the estimated road 
disturbance, is utilized to detect the current road level which the vehicle is driven on. The relationship between time- 
window and system response speed and estimation accuracy is studied and the optimal time-window length is determined. 
The influence of the variation of parameters in semi-active suspension on estimation accuracy is also studied. 

Two groups of tests based on a 1 /4 semi-active suspension test bench have been carried out to validate the real-time esti¬ 
mation effect in the actual driving process. In the first group of tests, the stability and adaptability of on-line estimation 
method of road profile were tested by changing the damping level and sprung mass respectively, and the results show that 
the method can well adapt to the variation of suspension parameters. In the second group of tests, a random sequence of five 
different levels of roads is utilized to verify the performance of the on-line estimation method by the research of road rough¬ 
ness classification. The test results show that the probability of correct estimation under all roads is more than 90%, and the 
overall accuracy of road classification reaches 98.5%. Furthermore, to validate the adaptability of on-line estimation method 
to the vehicle speed, the performance test under the same road level with different vehicle speed was carried out. The results 
indicate that this method can quickly adapt to the variation of vehicle speed, and the probability of correct classification 
reaches 96.25%. 

Although the on-line road profile estimation method designed in this paper performs wonderful performance in the road 
disturbance estimation and road roughness classification, it cannot truly reflect the comprehensive situation of the current 
road when a big elevation difference between left and right sides of the road (when a bump or a pit suddenly appears on one 
side of the road). Future research will consider the combination of road estimation and suspension control, design a semi¬ 
active suspension control method based on the road information perception, capable to coordinate the handling stability and 
ride comfort, and perform the test based on the 1/4 suspension test bench. 
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